Non-Markovianity assisted Steady State Entanglement 
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We analyze the steady state entanglement generated in a coherently coupled dimer system subject 
to dephasing noise as a function of the degree of Markovianity of the evolution. By keeping fixed the 
noise strength while varying the underlying dynamics, we demonstrate that non-Markovianity is an 
essential, quantifiable resource that may support the formation of steady state entanglement whereas 
purely Markovian dynamics governed by Lindblad master equations lead to separable steady states. 
This result illustrates possible mechanisms leading to long lived entanglement in purely decohering, 
possibly local, environments. We present a feasible experimental demonstration of this noise assisted 
phenomenon using a system of trapped ions. 



The generation and ultimate persistence of quantum 
entanglement is normally thought to be correlated with 
a high degree of system isolation, while the presence of a 
surrounding environment tends to decohere the quantum 
system thus driving it towards a classically correlated 
state. However, driven systems, generally out of equi- 
librium, have been shown to tend towards steady states 
where quantum correlations are nonvanishing [1]. A va- 
riety of entanglement preserving mechanisms have been 
put forward, showing how the presence of environmental 
noise can be instrumental in keeping the system entan- 
gled in the steady state [J-Q. In the framework of a 
Markovian dynamics, the presence of local pure dephas- 
ing is normally detrimental for entanglement preservation 
0. The situation though is very different when the de- 
phasing noise is non-Markovian. We show that the pres- 
ence of steady state entanglement can be unambiguously 
linked to an increasing degree of non-Markovianity, even 
when the environments are acting locally. This result 
is of particular interest in the light of identifying mecha- 
nisms that assist entanglement preservation in condensed 
matter and biomolecular systems, where non-Markovian 
dephasing is a dominant noise source 



The system - We consider a dimer system coherently 
coupled via an exchange interaction of strength J, so that 
the system Hamiltonian reads H s = tJ i cr / (T 7 
J{a^o~2 + o-^a^) (h — 1), and subject to the action 
of an environment that leaves the populations unaffected 
but tends to randomize the phases of superposition states 
(the so-called pure dephasing or transverse decoherence). 
We will model this situation by subjecting each element 
j = 1,2 of the dimer to the action of a localized mode, 
being the system-mode interaction governed by a Hamil- 
tonian of the form tL s -m = 9jO~iX a j + a ])> where aj(aj) 
denote the the operators of annhilitation (creation) of 
mode excitations. We will assume that the local modes 
are damped by a conventional Markovian bath so that 
the global time evolution of the dimer and the vibra- 
tional modes is described by a Markovian master equa- 
tion with a Liouvillian part accounting for the damping 
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FIG. 1: Beating pattern in the dimer population inversion for 
different values of / (See main text for details). A Markovian 
environment (/ S> 1) will simply wash out the discreteness 
of the vibrational modes while the presence of beating in the 
dimer signal is a signature of the persistence of a coherent 
interaction with a localized vibrational (damped) mode. 



of the localized modes at a certain rate n. Note, however, 
that tracing out the local mode leads to a density matrix 
for the dimer system that is in general no reproducible 
from a purely Markovian evolution for the dimer alone. 
We will show that the interplay between the entangling 
exchange interaction and the local dephasing will result 
in an entanglement dynamics that depends crucially on 
whether of not the action of the environment can indeed 
be described by merely subjecting the dimer to Marko- 
vian dephasing. In the limit where the coherent coupling 
g is smaller than the decay k we find that the effective 
dephasing rate j c g, which quantifies the strength of the 
system-environment coupling, is proportional to the ra- 
tio g 2 /n (see 



14j for the case of cavity QED). When 



2g <C k the dephasing is exponential and the dynamics 
of the system can be reproduced by a Lindblad master 
equation with the effective decay rate 7 c ff. For 2g < k, 
a regime which we denote as weakly non-Markovian, the 
definition of an effective dephasing rate 7 e fi is still reason- 
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able for times exceeding n^ 1 when the decay is well ap- 
proximated by an exponential. In the limit where g > k 
significant coherent oscillations can occur and the defini- 
tion of an effective dephasing rate becomes meaningless. 
In this work we will always be operating in the regime 
where 2g < k and will analyze the different dynamics 
that the dimer will undergo when subject to the same 
effective noise strength but experiences a different un- 
derlying dynamics. This will allow us to single out the 
precise influence of environmental time-correlations on 
the persistence/absence of stationary quantum correla- 
tions in the dimer system. 

A convenient parametrization- In order to link the 
observed dynamics to a parameter that quantifies in some 
form how much the environmental action departs from 
strict Markovianity, we will introduce an index / such 
that the coherent spin-mode coupling g and the mode 
damping rate k are given by, respectively: g = ^/fgo and 
k = fno. In this way, the noise strength g 2 / ' k — g^/no is 
kept fixed while varying / from values much larger than 1 
to much smaller than I leads to, respectively, Markovian 
and non-Markovian dynamics with, crucially, the same 
effective local dephasing rate 7 ff. In other words, the 
strength of the dephasing noise is kept fixed, but the un- 
derlying noisy dynamics is modified, making the coherent 
system-oscillator coupling to dominate over losses (small 
/ domain) or viceversa (limit of large /). Note, however, 
that while the decoherence rate is kept constant as the 
parameter f is varied, the population decay (which de- 
pends on g and k) does change, as shown in Fig. I. The 
monotonic change in Markovianity of the noise with the 
parameter f is shown in Fig. 2, where an explicit measure 
of non-Markovianity is considered. 

Analytical Results.- The effective Hamiltonian for 
the dimer-local damped modes system takes the form 

H = LOl<7^ cr-f + U2U2 <^2 + J( a t a 2 + °T a t) (1) 

+ (f2i — iK\)a\ai + (H2 — iK2)a\a2 

+gia[ z '(a 1 + a[) + g 2 (r£' (a 2 + a\), 

where uji and Qi (i = 1, 2) denote the site and the mode 
frequencies, respectively, and 

p = —iHp + ipH^ + 1n\a\pa\ + 2n2a2pa\- (2) 

Motivated by the situations frequently encountered in 
biomolecular complexes where external illumination is 
either weak and/or doubly excited states are strongly 
suppressed (Tsj . we will focus here on situations where 
the dimer dynamics is confined to the one excitation sec- 
tor so that the entire dynamics for the two sites takes 
place in the subspace spanned by the states |01), |10). 
By introducing the convenient delocalized states \u) = 
^(|01) + |10>) and \d) = -^(|01) - |10)) and the opera- 
tors 

<x p , m - -^=(a[ z) ±ai z) ), (3) 



a p , m = -/=(ai±a 2 ), (4) 

noting that o~ p \u) — a p \d) — while o~ m \u) — y/2\d) and 
o~m\d) = -\/2|w), in the subspace spanned by \u) and \d), 
the Hamiltonian part accounting for the system-mode 
interaction is given by 

H s - m = ^-a m (a p + a P ) + ^-a m {a m + a.J„), (5) 

where g± — g\ ± g 2 ■ Note that the restriction to the one 
excitation sectors removes the contribution from terms 
involving this operator does not couple to the 

states I it) and \d). In the special case where u\ = 102 = 
lj 1 f2i = O2 = O, Ki = K2 = k and 31 = 52 = g, the total 
Hamiltonian takes the form 

H = (u + J)\u)(u\ + (uj-J)\d){d\ (6) 
+(n - in)a\ n a m + V2g(\u)(d\ + \d)(u\)(a m + a* J, 

and the time evolution reads, 

p = -iHp + ipH^ + 2na m pa\ n . (7) 

An analogous derivation could be done for the case where 
both subsystems in the dimer couple to a global mode. 
In this case, H s - m = cr m g(a + a^), and therefore the 
dimer would remain decoupled if g\ — g2, while different 
couplings would yield to a system's dynamics analogous 
to that of the effective qubit discussed above 

with g = gi - g 2 . 

Figure 1 depicts the time evolution of the dimer's pop- 
ulation inversion of mer = 1 10) (10| - |01) (01| for different 
values of /. For this example the following parameters 
were considered: lj = 0, J = 1 (interdimer coupling), 
Q=2J, go = J and k = 20 J. With this, lcS = J/10. We 
truncated each oscillator after Fock layer 2 as the prob- 
ability to have a single excitation in either mode never 
exceeds 10% in the parameter regime that we consider 
here. Varying / in the range from 10" 3 to 100 allows 
us to move from a situation where g -C k to the do- 
main where g ~ k, so that the damped mode ranges 
from being Markovian to imprinting an element of non- 
Markovianity to the dynamics, as quantitatively exem- 
plified when evaluating the degree of non-Markovianity 
2? introduced in [18|. This measure provides a nec- 
essary and sufficient condition for a given evolution to 
depart from strict Markovianity by means of evaluating 
whether or not the associated dynamical map £( t + e ,t) is 
completely positive (CP) for any e. For that, it has to be 
(£(t+e,t) ® l) \®)(®\ > 0, where |$) is maximally entan- 
gled of our open system and some ancillary system [l9j . 
Using this condition, one can define a quantitative mea- 
sure to quantify how much does the evolution departs 
from strict Markovianity by evaluating X = J g(t)dt, 
where the function g quantifies how much does the norm 
of state (£( t+M ) ® 1) \&)(&\ differ from 1 in the limit of 



3 








qI , , , , 1 

20 40 60 80 100 120 140 

/ = e«/»/300 



FIG. 2: Using the non-Mar kovianity measure introduced in 
we can rigourously quantify the deviation from strict 
Markovianity of the dynamical map describing the time evo- 
lution of the dimer system. A finite value of the measure 
T>nm indicates that the dynamics cannot be accounted for 
in terms of a master equation of the Lindblad form. Here 
the values of / range from / m ;„ = 0.0035 to f ma x = 3.6554 
when the coordinate x varies between and 140. The value 
of T>nm decreases monotonically with / for the selected pa- 
rameter regime. For / sufficiently large, the dynamical map 
is divisible, T>nm is zero and the dimer's evolution is fully 
Markovian. 

e — > [18]. In Figure 2, we have plotted the normalized 
measure T>mm = yxi f° r a range of values of /, showing 
that the domain of large / leads to the dimer's dynamics 
to be fully Markovian. Limiting the harmonic oscillator 
to just 2 levels, we find for the steady state population 
of the state \d) the result 

4ff 2 + K 2 + (2J + ») 2 

In the limit of sufficiently small /, and in the regime 
where the local mode is quasi- resonant with the dimer 
eigenstates, so that il ~ J, one finds that the above ex- 
pression tends to unity. Hence, in this picture, the system 
is in weak interaction with a composite environment (lo- 
cal mode+Markovian reservoir) and relaxes towards its 
ground state, which in this case is the singlet state. When 
weakly coupled to a purely Markovian environment and 
experiencing dephasing at the same rate as before, the 
steady state of the dimer system tends towards the max- 
imally mixed state and no entanglement survives in the 
long term. The time evolution of the entanglement in the 
dimer system, as quantified by the logarithmic negativity 
of the state, is summarized in Figure 3, for an initial 
factorized condition where there is an excitation in one 
site only. Indeed, when / is large, the dynamics of the 
dimer can be reproduced by simulating a Markovian mas- 
ter equation with local dephasing rates 7 e ff. In this case, 
the steady state of the dimer is maximally mixed, a result 
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FIG. 3: Time evolution of the entanglement content of the 
dimer system when the index / varies from the range / 3> 1 
(Markovian evolution) to / < 1 (non-Markovian effects). 
When the evolution is strictly Markovian, the steady state 
is separable while in the presence of non-Markovianity, the 
dimer system contains quantum correlations that steadily be- 
come close to one e-bit for / sufficiently small. 

that can be proven using the results presented in section 
4 of [H, following theorem 5.2 of [171] • When / decreases 
but remaining in the weakly non-Markovian regime, the 
steady state approaches the singlet state, which in this 
case is the lowest energy eigenstate of the coherent evolu- 
tion. Provided that we keep on the weak non-Markovian 
regime (so that the condition m > g t is satisfied), alter- 
ing the symmetry of the problem and considering g± =/= §2 
does not reduce the overlap with the singlet below 90% 
up to ratios (72/91 = 2. At a finite temperature, the 
steady state entanglement is decreased but remains fi- 
nite. For an average photon number of n t h = 0.1, which 
corresponds to a temperature of about 100 mK for GHz 
frequencies, we obtain that the steady state entanglement 
reaches the value 0.85 for / = 0.01 (Finite T results are 
not shown in the Figures). 

Possible experimental realization - In this section we 
outline a possible implementation of the dimer in contact 
with non-Markovian environments that was discussed in 
the previous section. This realization employs proven el- 
ements of ion trap technology and offers, in principle, 
the possibility for control of all the relevant parameters 
in the Hamiltonian Eq. ©. Three basic dynamical el- 
ements have to be generated. The electronic degrees of 
freedom of trapped ions constitute the dimer and need 
to be coupled to each other, the motional degrees of free- 
dom play the role of the environment and need to be 
coupled to the electronic degrees of freedom via a disper- 
sive interaction and laser cooling needs to be employed 
to induce damping of the motional degrees of freedom 
which generates a Lorentzian spectral density. 

For the implementation of the direct coupling between 
the two qubits we make use of a S0rensen-M0lmer gate 
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2l[ l22j |. To this end one applies laser fields with two 
different frequencies so that the two-photon process cou- 
pling \gg) -f-> |ee) and \eg) -f-> |ge) is resonant, i.e., 
cjf + ^2 = 2w eg , while neither of the frequencies are 
resonant with single excitations of the ions. This can be 
achieved by choosing u>\ — u) eg — 6 and = ^eg + 8, 
where 5 <C ?7^_r, v. Here 77 denotes the Lamb-Dicke pa- 
rameter, flji the laser Rabi frequency and v is the fre- 
quency of the center-of-mass mode. If we prepare the 
two ions representing the dimer initially in a single exci- 
tation subspace we realize in this fashion the Hamiltonian 
-^dimcr = J c s(vt cr 2 + (J 2 (J i ) witn an effective exchange 
coupling J c ff ~ (flur/) 2 /(v — S) 21]. The added advan- 
tage of this scheme is the fact that it does not require 
the center-of-mass mode to be in its ground state thus 
increasing the robustness of the scheme. 
The excitation preserving coupling Hamiltonian a z (a + 
a^) between an ion and the motional degrees of free- 
dom can be achieved in several ways. One approach 

20] subjects the ion to a far off-resonant standing wave, 
which creates the state dependent potential V{z) — 
Vq cos 2 (kz+j)a z . Here k is the wave- vector of the stand- 
ing wave lasers. The operator z is readily expressed in 

terms of phonon operators, z = J2 n -^™^/5^~( a «+ a n) > 
where m is the ion mass, and M n is the amplitude of each 
vibrational mode n at the ion. Expanding this potential 
in the small parameter kz <C 1 we find that the leading 
order contribution is linear in z with the second order 
contribution canceling. Hence we obtain a Hamiltonian 

H sb = <J Z J2 n M n \j 2^""( a « + a ™)- Traveling wave fields 
can also be used to generate this type of interaction in a 
suitably transformed basis [20| . This scheme couples to 
all the modes of the ion crystal which may lead to correla- 
tions between the environments acting on separate parts 
of the dimer. We may select specific environment modes 
by choosing light fields that are slightly off-resonant to 
specific modes while being more strongly detuned from 
the remaining modes. 

Finally our dimer model assumes that the modes coupling 
to the constituents of the dimer are damped. Damping of 
motional degrees of freedom can be achieved, of course, 
by applying laser cooling to an auxiliary third ion that 
couples to all the motional modes of the system. To be 
effective, this auxiliary ion should be placed at the end 
of the ion chain to ensure coupling to all modes in the 
ion string. Combining these three elements yields the dy- 
namical Hamiltonian discussed in Eq. ([2]) which underlies 
the non-Markovianity driven steady state entanglement 
in a dimer system. 

Conclusions - We have demonstrated that, when 
subject to system-environment interactions of the same 
strength, the non-Markovian character of the noise can 
be the crucial property that leads to steady state en- 
tanglement where purely Markovian noise would result 
in the complete destruction of entanglement. A possi- 



ble experimental verification of this depahsing-assisted 
phenomenon in ion trap physics which employs only ex- 
perimentally demonstrated building blocks has been dis- 
cussed. We expect these studies to contribute towards 
the identification of the physical mechanisms that could 
underpin the persistence of stationary quantum corre- 
lations in very noisy environments occurring in natural 
conditions |9l-ll 31] . A key issue in this context would be 
the evaluation of this effect for realistic (possibly strongly 
non-Markovian), multicomponent, biological systems op- 
erating at physiological temperatures. Numerical meth- 
ods able to incorporate temperature dependence and the 
ability to simulate exactly structured spectral densities 
would be required to fully address this issue and initial 
steps towards this development have already been pre- 
sented [Hj|. 
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